Scaling and universality in turbulent convection 
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Anomalous correlation functions of the temperature field in two-dimensional turbulent convection 
are shown to be universal with respect to the choice of external sources. Moreover, they are equal to 
the anomalous correlations of the concentration field of a passive tracer advected by the convective 
fiow itself. The statistics of velocity differences is found to be universal, self-similar and close to 
Gaussian. These results point to the conclusion that temperature intermittency in two-dimensional 
turbulent convection may be traced back to the existence of statistically preserved structures, as it 
is in passive scalar turbulence. 

PACS numbers: 47.27.-i 



Heat and momentum transport in slightly heated flows 
are governed by the Boussinesq equations [1] 

dtv + V - Vv ^ -Vp - (3Tg + lyAv , ^ ' 

where T is the field of the temperature fluctuations, v 
is the velocity field, g is the gravitational acceleration, 
(3 is the thermal expansion coefficient and k, v are the 
molecular diffusivity and viscosity. The system is kept in 
a statistically stationary state by the external source of 
fluctuations /t. We will focus on the statistical proper- 
ties of temperature excursions at scales larger than the 
Bolgiano scale Lb - where buoyancy forces balance the in- 
ertial ones in the velocity dynamics - yet smaller than the 
forcing correlation length, L. In that range, temperature 
fluctuations cascade toward the small scales where they 
are eventually dissipated by thermal diffusivity. Dimen- 
sional arguments based on this phenomenological picture 
would lead to the Bolgiano-Obukhov scaling, in the range 
< r < L: 5^(r) = {[T{r,t) - T(0,t)]") - r-"/^, and 
S^{r) = ([(i;(r,t)-i;(0,t))-r]") ~ r^"/^ (see, e.g., Ref. [2] 
and references therein). Actually, due to the presence of 
structures of warm-rising or cold-descending fluid - the 
thermal plumes (see Fig. 1) - the statistics of temper- 
ature increments exhibits a nontrivial scale-dependence. 
Indeed, as shown in Figure 2, moments of temperature in- 
crements display a scaling behavior S'J(r) ~ r'»" charac- 
terized by exponents deviating from the dimensional ex- 
pectations. (Conversely, moments of velocity increments, 
~ ^ do '^ot show measurable deviations from 
dimensional scaling, i.e. Cn ~ 3n/5.) This anomalous 
scaling is a feature shared by a large class of turbulent 
systems: understanding the origin of this phenomenon 
from first principles is a major challenge of turbulence. 
With a significant exception. Indeed, for passive turbu- 
lent transport this problem has been recently solved. Let 
us briefly recall the main points, referring to Ref. 3 for 
a comprehensive review. We consider an idealized ex- 



periment of turbulent dispersion of a passive tracer, e.g. 
dye, in a convective flow. The equation that governs the 
dynamics of the concentration of tracer is 

dtC + v-VC^ kAC + fc , (2) 

while the velocity field v evolves according to Eq. (1). 
Although seemingly similar, the dynamics of tempera- 
ture and concentration fields are radically different: tem- 
perature is an active scalar, since it affects velocity via 




FIG. 1: Snapshot of the temperature and velocity fields. Dark 
areas identify cold regions. The Boussinesq equations (1) are 
solved in a two-dimensional doubly periodic domain, with 
1024^ collocation points. The Bolgiano scale Ib is compa- 
rable to the smallest resolved length-scale. Since in two di- 
mensions there is a net energy flux toward the large scales, 
a statistically steady state requires the momentum equation 
in (1) be supplemented by a friction term —av that drags 
energy from the gravest modes. As customary, diffusive and 
viscous terms are replaced by hyperdiffusive ( — kA^) and hy- 
perviscous (— i^A^) ones, to confine dissipative effects to the 
smallest scales. 
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the buoyancy forces, whereas concentration is a passive 
scalar. The concentration fluctuations show an anoma- 
lous scaling behavior S^{r) = ([C(r, i)-C(0, t)]") ~ r^S 
as well: the exponents differ from the dimensional 
expectation. Since S^{r) is a linear combination of 
various n-point correlation function of the concentra- 
tion field {C{xi,t) . . .C{xn,t)), the latter has to con- 
tain a contribution, denoted as Z^{xi, . . . , Xn), that car- 
ries the anomalous scale dependence. In mathematical 
terms, {Xxi, . . . , Xxn) = {xi, . . . ,Xn)- The 

main point is that the function Z^ is characterized by 
a special dynamical property that distinguishes it from 
a generic scaling function. Let us remind the passive 
scalar equation (2) can be written in the equivalent form 
-^C = fc, where ^ stands for the total derivative along 
the particle trajectories defined by the stochastic differ- 
ential equations dX = v{X,t)dt + V^dW{t), where 
W{t) is brownian motion. The remarkable result is that 
jj'{Z^)x = 0, where the total derivative is performed 
following n particles advected by the flow, and the aver- 
age is taken over the ensemble of all trajectories. In plain 
words, Z^ is statistically preserved by the flow [5, 6, 7]. 
In the specific context of a Gaussian, 5— correlated ve- 
locity field - the Kraichnan model - this is equivalent to 
say that the functions Z^ are zero-modes of the Fokker- 
Planck operator for n-particle diffusion [8, 9, 10]. An im- 
portant consequence of statistically preserved structures 
is that the passive scalar scaling exponents are universal 
with respect to the choice of the injection fc, since the 
latter does not enter the definition of Z^. 
We now turn our attention back to the temperature field. 
What we have learnt in the passive scalar case suggests 




FIG. 2: Scaling exponents of temperature, Cn , and velocity, 

Cn- The straight lines are the dimensional predictions, n/5 for 
temperature, 3ri/5 for velocity. Notice that at orders larger 
than n = 8 the temperature exponents saturate to a con- 
stant value ("^ ~ 0.8 [4]. The errorbars are estimated by the 
rms fluctuations of the logarithmic slope. To ensure the sta- 
tistical convergence of high-order moments we collected three 
hundreds of snapshots of the fields, spaced by half of the large- 
eddy turnover time L/(v'^}^^^. 



to investigate the effect of the external forcing fx on the 
scaling exponents. In Fig. 3 we show that the scaling 
exponents of temperature fluctuations are the same for 
two different choices of injection terms fx- Therefore, we 
conclude that the exponents C,^ are universal properties 
of two-dimensional Boussinesq convection. The univer- 
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FIG. 3: The moments of temperature differences, S'^{r), 
for n = 2, 4, 6, as a function of the separation r. Note 
the parallelism between curves of the same order n, within 
the scaling range. The two sets of curves are generated 
by two different kinds of injection mechanisms. In the first 
case (x), /t is a random Gaussian forcing, with correlation 
{fT{r,t)fT{r' ,t')) = Frdr - r'\)5(t - t'), where Ft decays 
with the characteristic scale L (approximately one-fourth of 
the box-size); in the second case the system is driven 
by the term fx = ■ u, that mimics the effect of a mean 
temperature gradient on the transport of temperature fluc- 
tuations. Wc consider only the isotropic contribution to the 
statistics, by averaging over all directions of the separation r. 
For orders equal or larger than n = 8 all exponents collapse 
- within errorbars - on the saturation value CS. — 0.8. The 
curves have been multiplied by appropriate numerical factors 
for viewing purposes. The equality of the scaling exponents 
for the two types of forcing has been checked by computing 
the logarithmic slope dlnS'X(r)/lnr (not shown). 



sality of scaling exponents suggests the possibility that 
a mechanism similar to that at work in passive scalar 
turbulence might be present in turbulent convection as 
well. To further pursue this line of thought, we notice 
that in the case of passive scalars, it is the whole func- 
tion Z^ to be universal with respect to forcing, not only 
its scaling exponent. It is thus of interest to look at the 
anomalous part of the temperature correlation function 
{T{xi,t) . . .T{xn,t)), to check whether it is universal. 
This measurement is unfortunately quite difficult for two 
reasons. First, the correlation function depends on 2n in- 
dependent coordinates; even if we exploit the statistical 
symmetries of this function - translational and scaling 
invariance - and we limit ourselves to the isotropic con- 
tribTition - by averaging over all configurations differing 
only by a rigid rotation - there will still be 2n — 4 degrees 
of freedom; in the most favorable case, n = 4 (for n = 2 
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FIG. 4: The functions cTT{r/R,9) = S2,2{R,r)/Sj{R) (left and center), and (7c{r/R,9) = S§^2{R, r) / (R) (right), in polar 

coordinates < r/R < 1 and < ^ < 2it, whore d = cos~^{R ■ f). The color is white where the function is zero, black where it 
is equal to unity. The function has a minimum at the origin, a{0, 9) = 0, and a maximum a{l, 0) = 1 at r = it. 



deviations from dimensional scaling are not detectable), 
the configuration space has four dimensions, which makes 
it quite untractable. Second, the anomalous part of the 
correlation function is hidden among several other con- 
tributions: it can be extracted only by taking proper lin- 
ear combinations, as, for example, in the case of S'^{r). 
To circumvent, at least partially, those problems, we fo- 
cus on a particular observable, S2 2{R,r) = {[T{R,t) — 
T(0,t)]2[T(r,t)-T(0,t)]2), which'is still anomalous, yet 
it has a nontrivial geometrical content. Since for r = R 
it reduces to the usual S4{R) we can write its functional 
dependence as S'|^2(-^''') = Sj{R)aT{r/R,0) where 6 is 
the angle between the directions of r and R. Since the 
scaling exponent of Sj{R) is universal, the bottom line 
is whether the "angular" part (Jt{t/R,6) is universal as 
well. In Fig. 4 we show a plot of the function a for the 
two different injection mechanisms (left and center). The 
similarity between the two pictures points to the conclu- 
sion that the anomalous part of the correlation function 
is again universal. 

This result leads us to conjecture that statistically pre- 
served structures Z'^{xi, . . . ^x^) might exist also for 
temperature: it is natural to define them by the property 
-^{Z'^)x = 0, as in the passive scalar case. Notice how- 
ever that, since temperature is an active scalar, this def- 
inition does not automatically ensure the universality of 
Z'^. Indeed, even if the forcing does not appear explicitly 
in the definition of Zj, the statistics of the trajectories 
X{t) in principle depends on fx, via the action of T on 
V. Therefore, should we accept the existence of statisti- 
cally preserved structures, the universality of anomalous 
temperature correlations requires to postulate that the 
whole statistics of v be universal as well. In Fig. 5 we 
show that this is indeed the case. 

Statistically preserved structures for temperature fluc- 
tuations entail another interesting consequence: since Z"^ 
is defined entirely in terms of the (universal) statistics of 



particle trajectories, and those are the same both for tem- 
perature and for concentration, we expect that Z^ = Z^ . 
In Fig. 6 we show that the scaling exponents of temper- 
ature and concentration are equal, as expected. As for 
the "angular part" of Z^, it is quite similar to that of 
Z^' (see Fig. 4, center and right). This is a further indi- 
rect evidence for the existence of statistically preserved 
structures for temperature statistics. 

In conclusion, the global picture of scaling and uni- 
versality in two-dimensional turbulent convection is as 
follows. Velocity statistics is strongly universal with re- 
spect to the external driving: probability density func- 
tions of velocity fluctuations are self-similar, and close 
to a Gaussian distribution, independently of the choice 
of fx- This is most likely a consequence of the observed 
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5rV/<(5rV)'>'" 

FIG. 5: Probability density function of longitudinal velocity 
increments 5rV = {v{r,t) — v{r, 0)) • r, rescaled to their stan- 
dard deviation ((5r«)^)^^^. We show two sets of data obtained 
by driving the system with random Gaussian forcing ( x ) , and 
by /t = 73 • w {+)■ Here r = 0.2, inside the scaling range. 
At different r the rescaled pdf's collapse onto each other, as 
expected for a self-similar statistics. The Gaussian density 
function is shown as a dotted line, for comparison. 
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FIG. 6: Local scaling exponents of temperature (x) and 
concentration (A) fluctuations, Cn''^ (r) = din Sn''^\r) /dlnr. 
Temperature and concentration are driven by independent 
Gaussian random forcings. 



tivc scalar turbulence, and depend on the specific form 
of the feedback of the scalar field on the velocity. For 
example, in three-dimensional turbulent convection the 
Navier-Stokcs equations arc characterized by a direct and 
intermittent energy cascade, the whole velocity statistics 
might then be nonuniversal, and a new type of universal- 
ity might emerge. 
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universal Gaussian behavior of the inverse energy cascade 
in two-dimensional Navier-Stokes turbulence [11, 12]. In- 
deed, velocity fluctuations in two-dimensional convection 
also arise from an inverse energy cascade which is driven 
now by buoyancy forces. At variance with the usual 
Navier-Stokcs inverse cascade, the energy injection now 
is not restricted to small scales. Indeed, the energy input 
rate e{r) = Pg ■ {v{r,t)T{0,t)) grows with the scale as 
e(r) ^ r^^^. This scale-dependent input induces the ob- 
served scaling S'^(r) ~ (e(r)r)"/^ ~ j,3n/5_ Temperature 
statistics shows anomalous scaling. This stems from the 
existence of statistically preserved structures, whose ex- 
istence explains the observed universality of anomalous 
temperature correlation fimctions and the equality be- 
tween temperature and concentration anomalies. 
Let us point out that the equivalence of the statistics 
of an active scalar, as temperature, to that of a pas- 
sive scalar, as concentration, depends crucially on the 
universality of the whole velocity statistics found here. 
That could however be a nongeneric phenomenon in ac- 
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